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Lecture 1 



It is great pleasure for me to visit China and to take part in this school. I would like to 
thank the Nankai Institute of Mathematics and Professor Ge Mo-Lin for hospitality. 

My lectures will be devoted to the Quantum Inverse Scattering Method (QISM). 
The QISM is a direction in the theory of quantum integrable systems which starts its 
history from the summer 1978 when 3 groups: in Leningrad, USSR (Faddeev et ai.), 
Fermilab, USA (Thacker, Creamer, Wilkinson) and in Freiburg, Germany (Honerkamp 
et aJ.), studying the quantum nonlinear Schrodinger equation, had found some striking 
and puzzling connections between the famous Bethe Ansatz and no less famous Classical 
Inverse Scattering Method (CISM). 

In other words, QISM has arisen as a result of a synthesis of two traditions in the 
theory of integrable systems which, up to 1978, were developing quite independently. 
The first tradition originates from the famous paper by Hans Bethe (1931) from which 
the Bethe Ansatz took its name. It was developed in the works of L.Hulthen, E.Lieb, 
C.N.Yang, C.P.Yang, R.J. Baxter and many other scholars devoted mainly to exactly 
soluble models of lattice statistical mechanics and of quantum mechanics. The second, 
more young, tradition originates from the paper of Gardner, Green, Kruskal and Miura 
on the KdV equation which gave rise to CISM (Lax, Ablowitz, Kaup, Newell, Zakharov, 
Shabat, Faddeev,... ). Still, the roots of this tradition can be traced even deeper into 
the mathematics of 19th century (Liouville, Jacobi, Kowalewsky,...). 

One should mention also two more directions which contributed considerably to 
QISM. These are the Factorizable -S-matrices theory due to A. B. Zamolodchikov and 
Al. B. Zamolodchikov and the various group-theoretical approaches (Adler, Kostant, 
Olshanetsky, Perelomov, Semenov-Tian-Shansky, Reyman et ai.). 

I am not going to present here the complete history of QISM. However, I would like 
to mention the main successes of QISM: exact quantization of the sine-Gordon equa- 
tion (Faddeev, Takhtajan) and calculation of correlators for various quantum integrable 
models (V. E. Korepin, N. M. Bogoliubov, A. G. Izergin, F. A. Smirnov). Witnessing 
the contemporary Quantum Groups boom I cannot help reminding you that the whole 
QG business has arosen as a by-product of QISM, being now quite independent disci- 
pline. Let me add also one more remark concerning the development of QISM. There 
is a general trend in QISM which appears to become stronger last years: shift of stress 
from applications and related analytical questions to algebraic structures underlying the 
integrability. 

In my lectures I am going to give an elementary introduction to QISM and to touch 
also some more special questions (Functional Bethe Ansatz). Since I would hke to 
concentrate on the mathematical methods involved rather than applications of QISM, 
in the center of our attention there will be the only problem: calculation of the spectrum 
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of integrals of motion. Various approaches will be illustrated on the sole example: XXX 
magnetic chain. 

To begin with, let me discuss briefly the concept itself of integrability. In the clas- 
sical mechanics there is well-known definition of integrability due to Liouville (1855). 
According to Liouville, the classical Hamiltonian finite-dimensional system is called in- 
tegrable if it possesses a set of independent integrals of motion commuting with respect 
to the Poisson bracket 

4} = 0. 

The total number of the integrals of motion (including Hamiltonian) should be the half 
of the dimension of the phase space. 

The Liouville's theorem provides also a way of constructing the action-angle variables 
for an integrable system in terms of the multi- variable curvilinear integrals. Being quite 
good for theoretical studies, this construction, unfortunately, does not help to perform 
an effective integration of concrete models. In practice, one needs usually to resort to 
more special techniques such as CISM or some algebraic methods. 

The situation with the definition of the quantum integrability is even worse. The 
first thing which comes to mind is to mimic the classical Liouville's definition, namely, 
to require the existence of N (= number of degrees of freedom) commuting operators 

[/,-,4] = o. 

Unfortunately, it seems to be very hard to develop this idea up to a satisfactory level 
of rigor. The main obstacle is the difficulty with correct definition of functional inde- 
pendence of integrals of motion in the quantum case. It is certain that it is necessary 
to restrict somehow the class of allowed functions but, as far as I know, at the moment 
there is no consistent theory of that kind. To sum up, I must confess that I don't know 
any good definition of the quantum complete integrability. 

Let me take during these lectures the most pragmatical point of view: I shall call 
a quantum system "integrable" if it is possible to calculate exactly some quantities of 
physical interest, such as the common spectrum of commuting quantum integrals of 
motion or some correlators (in these lectures I'll concentrate on the first problem: the 
spectrum). The word "exactly" needs, of course, some comment. 

An excursion into the history of classical mechanics and mathematical physics shows 
that the concept of "exact solubility" has been changed during the course of time. The 
main trend was the permanent extension of the class of functions used. First, elementary 
functions, then integrals of them, then solutions to certain second order differential 
equations, elliptic functions etc.. The most general concept of exact solubility elaborated 
in 19th century is, in my opinion, the concept of separation of variables, that is, reduction 
of a multidimensional problem to a series of one-dimensional ones. In my 3rd and 4th 
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lectures devoted to the "Functional Bethe Ansatz" I'll try to show you that this old 
concept works quite effectively in the framework of QISM. 

Now let us start the systematic introduction to QISM. The basic idea of QISM, as 
it has become clear the last years, is purely algebraic. Its roots can be traced up to 
the very dawn of quantum mechanics. I mean the early matrix formulation of quantum 
mechanics due to Heisenberg which led to a purely algebraic treatment of the quantum 
harmonic oscillator (Heisenberg, 1925) and the hydrogen atom (Pauli, 1926). The idea 
is to include the commutative algebra of the quantum integrals of motion {/„} into 
some bigger algebra A. The space of the quantum states of the system in question 
is considered then as a representation (usually, irreducible) of that bigger algebra A 
whose elements produce, roughly speaking, the transitions between the eigenstates of 
the quantum conserved quantities {In}- The common spectrum of {/„} can be found 
then by purely algebraic means. 

In case of the harmonic oscillator there is only one integral of motion, the number of 
particles operator N, and it can be included into the Heisenberg Lie algebra A generated 
by N and three extra generators: /i(central element), a, 

[h, a] = [h, a^] = [h, iV] = 

[a, a^] — h 

[N, a] = -a [N, a^] = 

The procedure of finding the spectrum of N using the creation/annihilation operators 
a, is well known. 

In case of the hydrogen atom (Coulomb problem) the components of the angular 
momentum vector and the so-called Laplace vector form an algebra whose factor over 
the relation Hamiltonian^ const is the so(4) Lie algebra. This remarkable fact allows to 
determine the spectrum of the Hamiltonian by purely algebraic methods. 

In the above examples the bigger algebra A was, up to minor reservations in the 
hydrogen atom finite-dimensional Lie algebra. The fundamental peculiarity of 

QISM consists in using a new class of algebras to describe the dynamical symmetry 
of quantum integrable systems. These algebras are neither finite-dimensional nor Lie 
algebras. 

The algebras used in QISM are described in terms of the generators Taf){u), a,(3 e 
{1 . . .d} which can be considered as the elements of the square d x d matrix T{u) 
depending on the continuous parameter u frequently called the spectral parameter. The 
associative algebra Tr is generated then by the quadratic relations 
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d 

/3i,/32=l 
d 

= E ^a2/32(^)^ai/3iM-R/3i/32,7i72(^ - ^) "^^^ ^ 
/3i,/32=l 

where Ra^a2.i3ii32{'^) some given function ("structure constants tensor"). Using the 
matrix notation 

T = T®\d T =i(l®T 
the former relation can be written down in the compact form 

R{u - v)T{u)T{v) = f{v)T{u)R{u - v) (1.1) 

The "structure constants" R are required to satisfy the consistency condition (the 
well known nowadays Yang-Baxter equation) 

d 

/3i,/32,/33=l 
d 

= E -Ra2a3,/32/33('^)-Rai/33,/3i73('" + '^^)-R/3i/32,7i72('") ''^ 

or, briefly, 

Ri2(u)Ri3(u + v)R23(v) = R23{v)Ri3(u + v)Ri2(u) (1.2) 

(the notation is obvious). 

It is easy to see that the matrix trace t{u) of T{u) 

d 

t{u) = tTT{u) = ^ T^a{u) 
a=l 

forms a commutative family of operators 

[t{u),t{v)] = \/u,v 

which can be thus considered as integrals of motion of some quantum intcgrable system. 
Generally speaking, there can be other independent integrals of motion but in the case 
(i = 2, to which I'll restrict my attention in these lectures, t{u) turns out to be the 
maximal commutative subalgebra of 7r. 

So, given a solution R{v) of the Yang-Baxter equation one can define the quadratic 
algebra Tr. Given a representation of the algebra Tr one obtains a quantum integrable 
system whose quantum space is the representation space of Tr and the commutative 
integrals of motion are t{u). The main problem of QISM is to find their common 
spectrum and, possibly, correlators of some physically interesting operators. 
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The main steps of QISM can be summarized as follows: 



1. Take an R matrix. 



2. Take a representation of Tr. 

3. Find spectrum of t{u). 

4. Find correlators. 

The first step imphes solving the Yang-Baxter equation (YBE). Many particular 
solutions has been found by trial-and-error method. As for for the general theory of 
YBE, the major contribution has been made by V. Drinfeld who gave an axiomatics 
of QISM based on the concept of HopJ algebras. In Drinfeld's axiomatics the steps 1 
and 2 of our scheme are intertwinned inseparably. In the base of his theory lies the 
concept of the quasitriangular Hopf algebra whose representations produce particular 
/^-matrices. Drinfeld has constructed also an important family of quasitriangular Hopf 
algebras called Yangians y[Q] and parametrized by a simple Lie algebra Q. 

I don't intend to go further into details of Drinfeld's theory. For these lectures 1 have 
chosen another approach supposing that a solution to the YBE is given from the very 
beginning. This standpoint lies closer to the original form of QISM as it appears in the 
works of Leningrad group and is more convenient for applications and for discussing the 
main topic of these lectures — the spectrum of t{u) (Step 3). 

From this very moment I'll consider only the case d — 2 that is T{u) being a 2 x 2 
matrix 

'^"^ [C{u) D{u)) 

Moreover, I'll restrict my attention to the simplest solution R{u) of the YBE, i?-matrix 
of the XXX-magnet. corresponding to the Yangian 3^[s/(2)]. The solution is expressed 
in terms of the permutation operator V in the tensor product (8) 



Vx ® y = y ® X Vx, y e 



and reads as 



R{u) ^u + r]V 



/ a 



\ 



b c 
c b 



a 
b 
c 



= U + T] 

= u 
= V 



[1.3) 



a ) 

in the natural basis {e\ Ci, ei (8) 62, 62 (8) ei, 62 (8 e-i\ in ® C^. 

Let us discuss now the second step of our scheme: finding representations of 7r for 
given R{u). The algebra 7r possesses an important property which is called comulti- 
plication. Let Ti{u) and T2{u) be two representations of 7r in the spaces Vi and V2 
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respectively. Then the matrix 

d 

T{u) = Ti{u)T2{u) Ta^{u) = J2 Tl,al3{u)T2fi^{u) 

is also a representation of 7r in the space Vi ® V2 called tensor product of representa- 
tions 71(1*) and T2{u). The possibility to multiply the representations of Tr provides 
immediately an opportunity to construct infinitely many representations given a set of 
particular representations. Usually such basic representations are chosen to have some 
simple dependence (say, polynomial) on the spectral parameter u. Traditionally, in 
QISM such elementary representations are called L-operators and their product 

T{u)^Ln{u)...L2{u)Li{u) 

resp. monodromy matrix. Every new representation of this type gives rise to a new 
quantum integrable system. Putting aside the question of completeness of such repre- 
sentations, I would like to remark that this family of representations serves perfectly all 
the models important for applications. 

Let us turn now to our basic example, the algebra 7r corresponding to the XXX R- 
matrix described above. This algebra possesses the remarkable central element (Casimir 
operator) called quantum determinant 

A(u)=q-detT(u) = tr^^'^ (u - (u + ^) 

= tn^f (tx + f)T(n-f)l^ 

= A{u-^)D{u+^)-C{u-^)B{u + ^) 
= A{u + ^)D{u-'j)-B{u+'j)C{u-^) 
= D{u+^)A{u-^)-C{u + ^)B{u-^) 

Note that the quantum determinant respects the comultiplication 



q-detTi{u)T2{u) = q-detTi{u)q-detT2{u) 



The simplest, one-dimensional representation of 7r is provided by a constant number 
matrix K satisfying the identity 

[R{u),K®K]^0 V-u. 

Actually, this condition is fulfilled for any matrix K due to the SL{2) symmetry of 
the R-matrix. Note that the quantum determinant of K coincides with its ordinary 
determinant 

q-deiK = det K 

We shall see in the next lecture that such representations describe the boundary condi- 
tions for the integrable chains. 
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The next in complexity goes the L-operator L{u) which is hnear in the spectral 
parameter u. It is constructed in terms of three operators Si, S2, 83 or S3, S± = Si±iS2 
belonging to the irreducible finite-dimensional (dim = 2/ + 1) representation of the Lie 
algebra s/(2) 

L{u) =u + riY,Sa(^a= ^^^^ ^1^^ 1 S± = Si±iS2 (1.5) 
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[Sa, Sp] — 



[S'3,S'±] — ±5"-!- 



2S- 
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Si + si + si = si + ^(5+5_ + = /(/ + 1) 

q-detL(M) = {u - Ir] - -){u + Ir] + -) 

Since the R-matrix R{u — v) in the quadratic relation ( |1 . 1|) determining the algebra 
7r depends only on the difference of the spectral parameters, it follows immediately that 
shift of the spectral parameter by a constant L{u) — > L{u — 6) is an automorphism of 
7r and provides thus a little bit richer family of representations. 

Now, the comultiplication property allows immediately to construct plenty of repre- 
sentations of 7r 

T{u) = KLm{u-6n)...L2{u-62)Li{u-6i)=(^^ ^ j (m) (1.6) 
The corresponding quantum determinant reads 

A{u) = detKY[ q-detL„(M - 5„) = det K Y[{u - - Lv - ^)(« - + Lv + 7^) 

n=l n=l ^ ^ 

To sum up, we have constructed a family of finite-dimensional (dim = n^=i(2/n + 1)) 
representations of Tr 



T[u\K G Mat(2, 2), iV G Z+, G Z/2}^^i, {5„ G C} 



n=\i 



parametrized by the matrix K, number iV of L-operators, spins /„ and shifts 5„. By the 
reasons which will be explained in the next lecture the corresponding quantum integrable 
system is called inhomogeneous XXX spin chain. 

As a matter of fact, the representations constructed turn out to be irreducible for 
almost all values of 5n ■ Moreover, the whole family contains all the irreducible finite- 
dimensional representations of the Yangian y[s/(2)], see prof. M. Jimbo's lectures on 
the present school. Let me remark that Yangian case corresponds to -ft^ = 1 only. 

In the next lecture we shall study these representations in more details. 
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Lecture 2 



In the last lecture we have constructed a huge family ( |1.6D of representations of 
the algebra Tr associated to the s/(2)-invariant R-matrix ( |1.3|) . What I am going to do 
now is to show that some of the representations constructed, namely, the representations 
differing by the order of factors Ln{u—6n) in the product (|1.6| ), are equivalent. Obviously, 
it is enough to prove the statement for the products Li{u)L2{u) and L2{u)Li{u). In 
fact, it follows from a general theorem of Drinfeld about the universal R-matrix claiming 
the equivalence Ti{u)T2{u) ~ T2{u)Ti{u) for any two representations Ti^2{u) of Yangian, 
but it seems me instructive to show the equivalence manifestly on the simplest example. 

However, let us consider first more simple representations K which are absent in 
Drinfeld's theory. It is obvious that the products K1K2 and K2K1 are not equivalent 
unless K^s commute. However, for the products of K and L situation is better. 

Theorem 2.1 IfdetK^O then KL{u) ^ L{u)K . 

Proof. It is necessary to find such an invertible operator K, in the representation 
space V of L{u) that 

KL{u) = lC-^L{u)KK. 
Let us use the s/(2)-invariance of the L-operator ( |1.5| ): 

[L(M),S" + ia,] = 0. (2.1) 

Since the matrix K is invertible it can be represented as an exponent 

^ 1 

a=l ^ 

Define now /C as 

3 

K. = ko exp ^ kaS". 

a=l 

From there follows immediately the identity 

[L{u),KIC] = 

proving the theorem. 

It follows from the theorem that one can transform equivalently any product of L's 
and invertible K's to the form (|1.6|) K being the product of i^'s in the same order. 
Let us return now to the products of L's. 

Theorem 2.2 Let two L-operators ( \1.3i ) represent the algebra in the spaces Vi and 
V2 and he characterised by the spins li^2 and shifts 61^2 respectively. Then an invertible 
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operator TZ in the space Vi intertwining the products Li{u — 6i)L2{u — 62) and 

L2{u - 52)Li{u - 5i) 

nU{u - 5i)L2{u - 82) = L2{u - 52)U{u - Si)n (2.2) 

and providing therefore the equivalence of representations 

Liiu - 6i)L2{u - ^2) ^ L2{u - 62)Li{u - 61) 

exists if and only if the difference 621 = 62 — 61 does not take one of the following values 

621 ^ ±{\l2 -li\ + IH Mh -h\+ 2)r/, . . . , ±{12 + h)v (2.3) 



Proof. Substituting the expression ( |1.5|) for L into equation (p7^) for TZ we obtain 
the equation 

niu -61 + 5f(T,)(M -62 + S^afs) = {u-62 + S^afs){u - 5i + 5fa,)7^ 

(here and further the summation over repeated indices is always supposed). Separating 
the terms containing u one notices that TZ should be S'L(2)-invariant 

[TZ, 5f + ^2"] =0 Va 

and depend in fact on the difference = 62 — 61 only. The remaining terms result in 
the equations 

7^[52l(Sf - ^2") + 2ie"^^5fS2^] = [521 (Sf - ^2") - 2ie'"^^S^S^]TZ Va (2.4) 
The S'L(2)-invariance of TZ allows to look for TZ in the form 

^= P.('^2i)P, (2.5) 

i=i;i-;2i 

where Pj are the projectors corresponding to the expansion of the tensor product of 
two finite-dimensional irreducible representations of SL{2) labelled by spins /i^2 into the 
sum of irreducible representations labelled by spin j. 

Using the last expansion ( |2.5| ) together with the equations ( |2.4| ) for TZ one obtains, 
after some calculation, the recurrence relation for the eigenvalues Pj (^21) of TZ 

,r . 621 + r/(j + 1) , . , 

= ^^^_^(^.^i) P.fe) (2.6) 

which determines TZ up to an insignificant scalar factor. It remains to notice that 
the necessary and sufficient condition for pj to be nonzero and therefore for TZ to be 
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invertible is the condition ( p^.3| ) for 61^2 and /i,2- The same condition, in fact, ensures 
the irreducibihty of the product Li{u — 6i)L2{u — 62) (for more details see Prof. Jimbo's 
lectures at the present school). In what follows the nondegeneracy condition (|2.3|) is 
always supposed to be fulfilled. 

Let us turn now to the quantum integrals of motion t{u) = tTT{u). The last theorem 
shows that the spectrum t{u) does not depend on order of L's in the product ( |1.6|) . 
However, the explicit expression of t{u) in terms of the spin operators Sn is not of 
course invariant under arbitrary permutation of Sn- 

For physical applications it is convenient to think of Sn as of spins of some atoms 
arranged as a one- dimensional closed chain (ring), so that it is natural to use the peri- 
odicity condition n + = n in all the formulas. The integrable models corresponding 
to the representation ( |l.(j| ) of the algebra 7^ would look more realistic if one could ex- 
tract from t{u) the local integrals of motion that is the quantities H^'^^ A; = 1, 2, 3, . . . 
expressible as the sums 

N 

Zlik) _ TT{k) 

~ -"n,n-l,...,n-fe+l 
n=l 

(the periodicity + 1 = 1 is supposed). The local densities should 
involve only k adjacent spins Sn, Sn-i, • • • , Sn~k+i 

[^S-i,...,n-fc+i> S^]=0 m>n or m<n + k-l 

An important case when such local integrals exist is that of the homogeneous spin 
chain corresponding to equal spins and shifts: ln = h 5„ = 5 = 

T{u) = KLn{u)Ln-i{u) . . . L2{u)Li{u) 

The homogeneous chain has the important property of translational invariance: the 
similarity transformation U{ )U^^ defined by the relations 

US^U~^ = ^^+1 USnU-' = /CiSi/Cr' (2.7) 

leaves t{u) invariant: Ut{u)U^^ = t{u). The invariance of t{u) follows directly from the 
cyclic invariance of the trace and from the theorem The transformation U{-)U^^ 
generalizes the ordinary translation for the periodic chain {K = 1) to the twisted periodic 
boundary condition specified by the matrix K. Note that, generally speaking, 7^ 1 
in contrast with the case K = 1. 

The translational invariance of t{u) suggests that the local integrals H^'''' should also 
be translationally invariant: 

^ -'^n,n-l,...,n-k+l^ ~ -"n+l,n,...,n-A:+2 
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The simplest one-point density is 



I leave the proof of the translational invariance [f/, = and commutativity [H^^\t{u)] 
for you as an exercise and pass to the problem of two-point Hamiltonian H^'^\ 

The following criterion proposed by Sutherland (1970) will help us in our search for 

Theorem 2.3 If for some translationally invariant two-point density Hl^^_i there exists 
a one-point translationally invariant density Qn{u) such that 

[H21, L2iu)L^{u)] = L2{u)Q^{u) - Q2iu)L,{u) (2.8) 

then H^"^^ commutes with t{u). 

Proof. The following calculation shows how the terms cancel in the commutator 

{u)Ln-l{u) ..] = 
. . .+ Ln+l{u)Ln{u)Qn-l{u) - (m)Q„ (m)L„_i (m) . . . 

. . . Ln+l{u)Qn{u)Ln-l{u) - Qn+l{u) Ln{u) Ln-l{u) . . . 

To finish the argument it remains to use the cyclic invariance of trace and the 
translational invariance of H^'^^ to show the cancellation of the boundary terms n = 1 
and n = N which I leave to you as an exercise. 

The density H^^^ can be described in terms of the operator TZ{S2i) constructed in 
the theorem Revising the proof of the theorem one observes that the same TZ{62i) 
intertwines operators Li{u) and L2{v) 

n{u - v)Li{u)L2{v) = L2{v)Li{u)n{u - v) (2.9) 

Using the formula ( p.6| ) it is easy to see that pj+i(0) = —pj{0) and hence 71(0) is 
proportional to the permutation operator V in V^V. Let us normalise TZ{u) multiplying 
it by appropriate scalar factor such that 71(0) = Vi2- Then the wanted expression for 
H'21 is given by H^^ = 'R{0)Vi2 where TZ stands for the derivative dlZ/du. 

To prove the statement let us differentiate ( |2.9| ) with respect to u aX. u = v. Multi- 
plying the result 



H2iVi2Li{u)L2{u) + Vi2Li{u)L2{u) = L2{u)Li{u)Vi2 + L2{u)Li{u)H2iV: 



12 



from the right hand side by V12 and using the obvious identities V^V^ = 1, V^LiVu 
L2, Pi2-Z^2'^i2 = Li we arrive at the Sutherland equation ( p.8|) for Qniu) = Lniu). 
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Using the formula ( |2.(j| ) for the eigenvalues of TZ one can derive the eigenvalue ex- 
pansion for 

h+h 3 

H21 = const + coeff x ^ -^j' X! Z 

j=\h-l2\ k = i ^ 

One can show that the higher local integrals H^^^ etc. also exist for the homogeneous 
XXX spin chain. The most popular method to produce them belongs to Baxter (1972) 
and Liischer (1976) and gives an expression for H^^^ in terms of higher derivatives of 
TZ{u) at M = 0. I'll not go into details of the calculation which are well described in the 
literature and simply reproduce the result. 

Let us construct first the quantities t^''\u) which are easier to describe by words 
then to write a compact formula. Let t{u) = t{u + |). To construct t^''\u) take the 
product t{u)t{u + 1]) . . . t{u + {k — 1)77). Then subtract all possible terms obtained from 
it after replacing two adjacent factors t{u + kri)t{u + {k+ 1)1]) by A{u + {k + ^)ri) where 
A(m) = A(m + |) is obtained from the quantum determinant ( p..4|) by the same shift as 
t(u) from t{u). Then add all terms obtained from the original product after replacing 
in the same way two pairs of t's by A's and continue changing sign at each step until 
the t's will be exhausted. The result is t^''^{u + (A; — 1)|). For example 

t^^\u) = t{u) 
t(2)(n+|) = t(n)t> + r^)-A(^ + |) 

t^^) {u + T]) = t{u)t{u + 7])t{u + 2t]) - A{u + |)t(n + 2r/) - t{u)A{u + y ) 
(u) = t{u)t{u + ri)t{u + 2ri)t{u + 3r/) - A(m + |)t(M + 2r])t{u + 3r/) 
-t{u)A{u + y )t(M + 3r/) - t{u)t{u + ri)A{u + y ) 
+A(.. + |)A(^ + ^) 

and so on. 

It turns out that if I is the spin of the homogeneous XXX spin chain then t^'^^\u) taken 
at li = is proportional to the generalised translation operator U defined previously 
( |2.7|) . Finally, the local integrals of motion H^'^'^ are obtained as the coefficients of the 
power series 

00 k 

t{u) = ln[/-4(2')(«) = E TjH^'^'^ (2.10) 

k=i ^■ 

In conclusion, I want to show you a simple method for generating higher integrals 
which is not so widely known. The method is based on the concept of master symmetry 
which can be defined as an operator B producing one integral of motion from another 

i/^'^)] = (2.11) 
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Tetelman (1982) has shown that for the integrable spin chains the master symmetry 
B can be taken as the following sum ( "boost" operator) 



(2) 



n-1 



Strictly speaking, the above formula makes sense only for the infinite chain since 
the linearly growing coefficient n contradicts obviously the periodicity of the spin chain. 
Nevertheless, in the periodic case one can use it quite formally in commutators which 
"differentiate" the linear dependence on n and yield periodic expressions. 

The possibility to use ;B as a master symmetry is based on the obvious identity 

[B,U'']=nU''H^^^ (2.12) 

and the identity 

[B,t{u)]=i{u) (2.13) 



which is proven in the same way as the theorem |2.3| with the only difference that due 
to the factor n in the definition of B the terms do not cancel completely but add up to 

i{u) 

[... + {n + l)ifi5i,n + ^^S-i + •••,••• L„+i(u)L„(u)L„_i(u) . . .] = 
. . . + nLn+i{u)Ln{u)Ln-i{u) - nL„,+i{u) L„,{u) Ln-i{u) + ... 
. . . + (n + 1)L„+i(m)L„(m)L„_i(m) - (n + 1)L„+i(m)L„(m)L„_i(m) + . . . 

Using the identities (|2.12|) and ( p.l3| ) one obtains the commutation relation 

[B,r{u)] =t{u)-H^^^ 

between B and the generating function ( |2.1(]| ) of H^^\ Expanding the last equality in 
powers of u one arrives at the wanted relation ( |2.11| ). 



Lecture 3 

Now we are in a position to attack the main problem of QISM — the spectral 
analysis of the commuting family t{u), see Step 3 of the general scheme given in the 
Lecture 1. Traditionally, various methods of solving the problem have the common 
name: Bethe Ansatz (BA), with the corresponding adjective: coordinate BA, algebraic 
BA, analytic BA, functional BA, nested BA etc.. In the last two lectures I'll touch 
briefiy the algebraic BA (ABA) described at length in Prof. L. D. Faddeev's Nankai 
lectures (1987) and concentrate on the new variant of BA — the functional BA. 
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The ABA works for the highest vector representations T{u) of the Yangian 3^[s/(2)] 
that is the representations possessing the heighest (vacuum) vector |0) which is annihi- 
lated by the off-diagonal element C{u) of the matrix T{u) and which is an eigenvector 
of its diagonal elements A{u) and D{u) 

C(u)\0)-0 ^("^1°^ = Wu 

It is easy to show that the eigenvalues A±(m=i= |) factorize the quantum determinant 
A{u) of T(u) 

A+(«-|)A-(« + |) = AH 
Let us look now for the eigenvectors of t{u) in the form 

M 

\Vi, V2,..., Vm) = n B{Vm) |0) 
m=l 

Then it is possible to show that the eigenvalue problem 

t{u) \vi,V2, Vm) = r{u) \vi,V2, vm) 
is equivalent to the set of Bethe equations 

— T = 11 — m = l,...,M (3.1) 

A_(fm) ^:^iVm-Vn + r] 

the corresponding eigenvalue r(n) of t{u) being 

r{u) = A_iu) n """^"'"^ + A+(«) n """^"^"^ (3.2) 



Note that one can obtain the Bethe equations (|3.1|) from the formula (|3.2|) for t{u) 
by taking residue at u = Vm and using the smoothness of the polynomial t{u). The 
last equation ( p.2|) is in turn equivalent to the linear finite-difference spectral problem 
(Baxter's equation) 

t{u)Q{u) = A+(m)Q(m + ri) + A_{u)Q{u - r]) (3.3) 



M 

m=l 



for the polynomial Q{u) determined by its zeroes {vm} 

M 

Qi^) = n ~ ^rn) 
m=l 

We shall return to the Baxter's equation when speaking about the Functional Bethe 
Ansatz. 
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Let us apply the ABA construction to our favorite example — the inhomogeneous 
XXX magnetic chain. It is a simple excersize to show that the vacuum vector exists if 
and only if the matrix K determining the boundary condition is triangular {K21 = 0). 
For the sake of simplicity we shall assume K to be diagonal 



In that case the vacuum |0) is the common heighest vector of the local representations 



The corresponding eigenvalues of A{u) and D{u) are polynomials of degree N 



The ABA proved to be quite powerful method which helped to solve such models 
as sine- Gordon, XXX and XXZ magnets and others. However, being restricted to the 
heighest vector representations, it fails for representations which do not satisfy this con- 
dition. In particular, such interesting integrable models as sinh-Gordon, Toda chain, 
quantum tops fall beyond the reach of ABA. One should mention also the so-called 
completeness problem that is the question if the Bethe eigenstates \vi,V2, ■ ■ ■ ,vm) are 
complete. There are simple examples when this is not the fact and some special inves- 
tigation is needed in order to gain the missing eigenstates. 

In the rest of this lecture and in the last one I'll describe an alternative method, the 
Functional Bethe Ansatz (FBA), which is free of the restrictions inherent in ABA. 

The Functional Bethe Ansatz was born from a fortunate marriage of two ideas. The 
first one is the central idea of QISM that is using quadratic R matrix algebras (Yangians 
etc.) to describe the dynamic symmetry of integrable systems. The second one is the 
very old idea of separation of variables. So, before I'll proceed to FBA itself, let me 
remind you some elementary facts concerning the separation of variables. 

Briefly, the separation of variables is a reduction of a multidimensional spectral prob- 
lem to a system of one-dimensional (multiparameter) spectral problems. The simplest 
example of separation of variables is provided by the Cartesian coordinates on the plane 

R2 3 (xi,X2). 

Let us introduce two commuting differential operators of second order 




Sn of SU{2) 



S-\0)^0 SM^-ln\0) yne{l,...,N} 



N 



n=l 



dxi 



dxl 
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which can be thought of as the integrals of motion of some integrable quantum system. 
It follows from the commutativity that the operators Di possess common eigenfunctions 

It is well known that the eigenfunctions factorize 

$(a:i,X2) = Xi(xi)X2(x2) 

the original two-dimensional spectral problem being thus reduced to two independent 
one-dimensional spectral problems 



X'( — \iXi 
X2 = A2X2 



Less trivial example is presented by the parabolic coordinates (2/1,1/2) 

- - yl 

- 2 

= ym 

The corresponding commuting second-order differential operators are 



^ ^ 92 I 92 ^ 1 / 32 I 92 

'dxl fef 3/2+2/2 \dyl dy^ 

^ ^"^1 dx'^ ~^ '^'^'^ dxidx2 ~^ dxi y'f+y'i (^2 Qy'j Vl Qy'^ 

As previously, their common eigenfunctions 

A<l> = A$ D^ = fi^ 

factorize 

Hyi,y2) = Yi{yi)Y2{y2) 
The corresponding separated equations are 

r YI'-iXy' + fi)Y, = 
\ - (Ay2 - ^)Y2 = 

Note that, in contrast with the previous example, the spectral parameters A and fi 
cannot be decoupled, and both equations should be solved together. This situation of 
the so called Multiparameter spectral problem is typical for the separation of variables 
in the general case. Note that the Baxter's equation ( |3.3| ) arisen in ABA is another 
example of (finite-difference) multiparameter spectral problem. 

In the above examples the separation of variables was produced by a pure change of 
coordinates. In the most general case, however, the separation should be obtained via 
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generic canonical transformation involving not only coordinates but also momenta. The 
simplest example of this kind is the Goryachev-Chaplygin top. 

The quantum GC top is described in terms of the operators Ja 
belonging to the e(3) Lie algebra 

[Xa,Xi3] = 

The values of the Casimir operators are assumed to be fixed 

3 3 
a=l a=l 

Note that the second constraint is crucial for integrabihty. 
The Hamiltonian of the model is 

H=\{Jl + Jl + AJl) - bx, = + 3Ji) - bx, (3.4) 

where 6 is a parameter (magnitude of the external field) . 
The second integral of motion commuting with H is 

G = 2 J3( J2 - J| + ^) + 6(x3 Jl + J1X3) (3.5) 

(note the quantum correction |). 

The separated variables {u,, U2) in the classical case were found by Chaplygin himself 

Ui = Js + VJ^ 
U2 = J-i - 

Their quantum generalization has been found few years ago by Komarov. They are 
defined as the "operator roots" of the quadratic polynomial with commuting coefficients 

- 2J3M - ( J2 - Jl + ^) 

(note the quantum correction |). The common spectrum of (1*1,1*2) turns out to form 
an equidistant square lattice 

1 1 

spec(Mi, U2) = {1 1: - + 2ni, -1 ± - - 2n2} 

The common eigenfunctions of H and G factorize in the li-representation 

$(M1,M2) = ^{Ui)ip{u2) 

The corresponding one-dimensional separated equation is the same for u, and U2 

t{u)lp{u) = d{u + l)ip{u + 2) + d{u - l)^p{u - 2) 
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where 



3 



2{H+-)u-G 



d{u) 




The example of GC top shows how nontrivial can be the task of finding the separation 
of variables. There is no general recipe guaranteeing success for any given integrable 
model. Fortunately, for the integrable systems subject to QISM such a recipe, to which 
I have given the name "Functional Bethe ansatz (FBA)", exists, at least in the case of 
the Yangian y[sZ(2)]. The success of FBA lies in using the algebraic machinery of QISM 
which provides a powerful tool for finding the separated variables. 

So, I am starting the detailed exposition of FBA. Let me remind you that throughout 
these lectures we consider the algebra Tr associated to the R matrix of the XXX model. 
Suppose that the 2x2 matrix T{u) is a polynomial in the spectral parameter u 



and defines a representation of 7r in a finite-dimensional space W . 

The recipe which constitutes the heart of FBA can be expressed in few words. Take 
the polynomial B[u) having commuting operator coefficients (the commutativity 



the operator zeroes x„ of B{u) provide the wanted separated variables. The rest of my 
lectures will be devoted to the deciphering and adjustment of this obscure remark. 

The first problem which arises immediately is to give the precise meaning to the roots 
of the operator polynomial B(u). In order to simplify the task and to avoid troubles 
with degenerate cases I'll impose few conditions on the representation T[u). The first 
condition is always fulfilled for irreducible representations of Tr. 

Condition 1 The senior coefficient T/v is a number matrix. The quantum determinant 
A(m) ofT{u) is a number-valued function. 

The second one is a nondegeneracy condition. It ensures that B{u) and A(m) are 
polynomials of the maximal degree. 

Condition 2 




[B{u),B{v)] = Q Wu,v 



is a simple consequence of the basic quadratic relation ( |1.1| ) for the T matrix). Then 



Bn 7^0 



det Tn ^0 
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Though it is not yet clear what are the operator roots of B{u) their symmetric 
polynomials bn are easy to define 

bn = {-ITBm-u/Bn n = 1, 2, . . . , TV 

Using bn one can express B{u) as 

B{u) - 5jv(m^ - + hu''-'' ) 

It follows from the commutativity of B{u) that the operators {bn}n=i ^^^^ commute. 
Let B be their common spectrum. 

B = spec{6„}^=i 

The following condition of semisimplicity is the most restrictive one. Fortunately, it 
holds for many interesting models. 

Condition 3 The operators {bn]n=i have a complete set of common eigenfunctions and 
to every point h — {bi, . . . , bpf) G B C there corresponds only one eigenf unction. 

It follows from the last condition that the representation space W is isomorphic to 
the space FunB of functions on the finite set B C C^. Of course, the isomorphism 
is not unique. It is defined up to multiplication by some nonzero function on B. So, 
let us suppose that some isomorphism is fixed and the operators 6„ are realized as the 
operators of multiplication by the corresponding coordinates 6„ in 

b= e /eFunB 

M(b)=6n/(b) 

In what follows in such cases I will not make difference between the operators like 6„ 
and the corresponding functions like bn and always will omit hats over operators (see 
later operators 

Since 6„ are the symmetric polynomials of the roots of B{u) which are yet to be 
defined we are led to consider the mapping 

e : ^ : X ^ b 

given by the formula 

where s„(x) is the elementary symmetric polynomial of degree n e {1, . . . , iV} of N 
variables {xn\n=i- 

Sl(x) = X1+X2H \-xn 

S2(x) = X1X2 H 

Siv(x) = XiX2---Xn 
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Note that the pre- image ^(B) is symmetric under permutations of the coordinates 

{^n}n=l- 

And finally, the very last condition. 

Condition 4 Pre-image X = 6~^(B) contains no multiple points that is each b e B 
has exactly N\ pre-images. 

In principle this condition can be omitted but in that case one is urged to use rather 
sophisticated language of jets in order to work with the multiple points. 

It follows from the last condition that the mapping G induces isomorphism between 
W = FunB and the space SymFunX. Now it is easy to define the operator roots x„ of 
B{u) as multiplication operators in the extended representation space 

W = FunX dW^ SymFunX 
It is obvious that spec{xn}n=i — X and that 

iV 

B{u) = Sat n - Xn) (3.6) 

n=l 

For the last equality to be correct its right hand side should be restricted to the 
space W since B{u) is originally is defined only on W. However, the same equality can 
be considered as definition of the natural extension of B{u) from W to W. 

Working with W one should always keep in mind that this is a "nonphysical space" 
and all the final results should use only the original space W — FunB=Sym FunX. 

The operators x„ being defined the next problem is to calculate the expression for the 
commuting integrals of motion t{u) — trT{u) — A{u)+D{u) in the x- representation and 
to observe the resulting separation of variables. This task will be solved in several steps. 
First, let us introduce the "momenta" conjugated to the "coordinates" Consider 
the operators 

N 

=E<^n=[^HU.„ (3.7) 
p=l 

N 

X: = Y:<Dn^[Diu)U.^ (3.8) 
p=i 

which are obtained from the polynomials A(u) and D{u) as their values at the "points" 
u — Xn- The substitution of operator values for u will be defined correctly if one 
prescribes some rule for operator ordering. Here the ordering of x's to the left is chosen. 
I shall call it "substitution from the left" . Note that the operators X^ act from W to 
W. 

It is convenient to introduce the shift operators — er^^^/^^n 

: ^ : {xi, . . . ,Xn, . . . , xn) ^ (xi, . . . , a:„ ± r^, . . . , xn) 

It turns out that the operators X^ have nice commutation relations with Xn- 
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Theorem 3.1 

X^Xn = (Xn ± ri5mn)X^ Vm, 71 (3.9) 

Proof. Consider the foUowing commutation relation between A[u) and B{u) which 
can be extracted from the basic commutation relation (|1 . 1| ) 

{u - v)A{v)B{u) + r]B{v)A{u) = {u - v + r])B{u)A{v) 

Now substitute Xn for v from the left. The second term in the left hand side will 
vanish since B{xn) = 0. 

(U - Xn)X~B{u) = {U-Xn + 7])B{u)X~ 

After substituting the expansion ( p.6| ) for B{u), cancelling the factors B^ and {u — 
Xn), and expanding both sides in powers of u we obtain the relation 

X-s(x) = s(E-x)X- (3.10) 

first for any elementary symmetric and hence for any symmetric polynomial s(x). Note 
that, since X is a finite set, s(x) can be in fact any symmetric function on X. 

Quite similarly, the commutation relation between X^ and Xn is obtained from the 
identity 

{u - v)D{u)B{v) + r]D{u)B{v) = i^u - v + t])B{v)D{u) 

The next natural step would be establishing the commutation relations between 
X^'s. However, it cannot be done immediately because of the fact that the operators 
X^ act by definition from W to W . So, we need first to extend them from WioWas, 
we did it with B{u). To this end, consider the constant function on X 

a;(x) = 1 Vx e X 

Obviously, u is symmetric under permutations and thus belongs to SymFunX = W. 
Consider now the functions on X 

A±(x) = [X±.;](x) VxGX 

defined as the images of uj for the operators X^. 

Let us show now that due to the commutation relations ( p.9|) between X^ and x„ 
the functions A,^ on X determine uniquely the action of X^ on any vector s of W. 
Note that any vector s ^ W which is identified due to the isomorphism W = FunX 
with some symmetric function s{xi, . . . , x^) can be created from the cyclic vector u by 
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the operator s = s{xi, . . . ,xn)- Then, applying to s and using the identity ( |3.10| ) 
for the function s(x) we obtain 

[X±s](x) = [X±5^](x) = s(E±x)[X±a;](x) = 5(E;^x)A±(x) 

Now we can consider the last equality as the definition of action of X^ on arbitrary, 
not necessarily symmetric, function s G FunX = W. Obviously the commutation 
relations (|3.9| ) with x's will still be valid for X^'s thus extended on W. Consequent ely, 
the relation ( p.lOj ) is now valid for any (not necessarily symmetric) function s(x). In 
other words, the extended operators X^ can be expressed in terms of multiplication and 
shift operators 

Now we can calculate the commutation relations for extended X^'s. 
Theorem 3.2 

Vm, n 

\/m,n m ^ n (3.12) 
A(a;„ ± |) Vn 

Proof. Let us prove first the commutativity of X~'s. The assertion being obvious 
for m = n, it is enough to consider the case m = 1, n = 2. Take the product A{u)A{v) 
and substitute in it m = Xi and v = X2 

[A{u)A{v)]u=.. = ^ x^'xTAn.An, = ... 
Then, using the commutativity of x's and the definition of Xf we obtain 

. . . = ^ ^ X2^ X^^ A^j^A^^ = ^ ^ X^ An = . . . 

nin2 n 

Finally, using commutativity of X2 and Xf and the definition of X^ we arrive at 
Xf X2" 

. . . = ^ X^ X2 A„ = X^ X2 

n 

In the same way, starting from A{v)A{u) one obtains X2"Xf . Since A{u) and A{v) 
commute due to ( p..l[ ) so do Xf and X2', the first assertion of the theorem being thus 
proven. Quite analogously, the commutativity of X"'''s follows from the commutativity 
of D{u). Similarly, commutativity of X+ and X~ for m 7^ n is derived from the identity 

{u - v)D{u)A{v) + r]B{u)C{v) = {u - v)A{v)D{u) + r]B{v)C{u) 



n n 
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Consider now one of equalities ( |1.4| ) for the quantum determinant 

A{u - |) = A{u)D{u - r^) - B{u)C{u - t]) 

Substituting u = Xn from the left leaves only the first term in the right hand side. 
Proceeding in the same way as above we obtain 



A(xn - i) = [A{u)D{u - r])]u=.^ = E <(a;„ - V^ApD^ = E(x„ - vYxPApBg 
1 



Now it remains to use the commutation relation ( |3.9| ) between X„ and x„ to arrive 
at the wanted result 

... = J2x-xiD, = x-x:i 

1 

The analogous equality for X+X^ is obtained in the same way starting from the 
identity 

A(m + |) = D[u)A[u + r^) - B[u)C[u + i]) 

I propose as a home exercise to look once more the proof of the theorem through 
and to trace where the extended definition of X^'s was used and where the original one 



Let us write down the whole set of commutation relations between x's and X^'s 

[x„,x„] = Vm,n 

X^Xn = {Xn±r]Smn)X^ Wm,n 

[X±,X±] =0 Vm,n 

[X+,X^] = Vm, n m ^ n 

X±X^ = A(x„±f) Vn 

It is quite natural to inquire about the representation theory for the associative 
algebra Af^ defined by the generators Xj^}^^^ and the above relations, and labelled 
by the polynomial A(m) of degree 2X. 

In the case of finite-dimensional representations such that the spectrum X of the 
operators is simple and has no multiple points, as throughout my lecture due 

to the Conditions 1-4, the problem of constructing a representation of A'a is equivalent 
to that of finding the functions {A^}^^^ on X satisfying certain relations which follow 



from the relations (|3.12|) between X^'s 



A^(x)A^(E±x) = A±(x)A±(E„±x) Vm,n Vx G X 

A+(x)A;(E+x) = A;(x)A+(E~x) Vm,n m + n Vx G X (3.13) 

A^(x)AJ(K„±x) = A(x„±|) Vn Vx G X 
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Strictly speaking, the above relations are not defined when the shift moves the 
point X out of the set X. It means that the functions {A^}^^;^ must vanish on the 
boundary X^ of the set X with respect to the shift 

A±(x) = VxeX± 

(3.14) 

x± = {xex|£;±xe(c\x)} 

As a consequence, the function A(m) must vanish on the set 

u((x:+|)u(x„--|)) 

Irreducibility and finite-dimensionality must put severe restrictions on the set X and 
the functions A(m), AJ. It seems highly probable that a kind of Stone-von Neumann 
theorem should exist for the algebra Aa that is there should exist essentially unique 
irreducible representation for every allowed A(u). 

Let us consider a simple example of irreducible representation for Aa- Let numbers 
{y^n}n=i be such that (A+ — A~)/?7 = 2Z„ is a positive integer Vn. Let the sets A„ be the 
equidistant (step 77) strings connecting 

An = {A-, A- + 77, . . . , A+ - 77, A+} |A„| = 2Z„ + 1 (3.15) 

Let us suppose, in addition, that the sets A„ do not intersect 

m^n =^ A^nA„ = (3.16) 

Define now X as the parallelepiped 

X = Ai X A2 X • • • X Ajv |X 

and the functions A^, A(m) as 

A±(x) = A±(x„) 

JV 

A±(«)-e±n(«-^n) e±^o 

n=l 

A(«) = A_(« + |)A+(«-|) (3.18) 
where ^± are arbitrary nonzero numbers. 

Theorem 3.3 The functions A^ define an irreducible representation of the algebra Aa 
in the space FunX. 



= n(2/n + l) (3.17) 

n=l 
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Proof. One can verify easily that the functions satisfy the conditions ( ^.131 ) 



( p.l4|) and hence define a representation of X/\ in FunX. Its irreducibihty follows from 



the standard argument. Suppose V C FunX is an invariant subspace. Being invariant, 
in particular, under the commutative subalgebra generated by {xn}n=i the space V 
can consist only of the functions vanishing on some subset of X. However, such a 
subspace cannot be invariant under the operators because of the condition (|3.16|) 
which ensures that the functions A^(x) have no zeroes on X other than in X^. 

Theorem 3.4 Let functions define a finite- dimensional representation of with 
the same A(u) and X as for the above described standard irreducible representation. 
Then the two representations are equivalent. 

Proof. We are going to show that the equivalence is provided by the multiplication 
operator p : FunX FunX : /(x) p(x)/(x) where p(x) is a function having no 
zeroes on X. The equivalence of operators leads to the equation 

p(E±x)A±(x)=p(x)A±(x) 

which can be considered as a set of recurrence relations for p(x). Note that the recurrence 
relations are compatible because of the conditions ( ^.131 ) for A^ and A,^. The function 
p(x) is thus defined uniquely on X up to an unsignificant coefficient. It remains to 
show that p(x) 7^ on X which follows from the fact that the functions A^(x) have no 
zeroes on X other than on Xj. The last assertion follows in turn from the third of the 
equalities ( |3.13D , definition ( p. 181 ) of A(m) and the condition 



It is natural to ask whether the above constructed sample representations exhaust all 
the irreducible finite-dimensional representations of A'a. I don't know the answer and 
hope that somebody of the audience will be interested in this problem. Fortunately, the 
results already obtained are enough to help solution of our main problem — spectral 
analysis of r(u). 

Lecture 4 

Let us apply now the results of the previous lecture concerning the representation 
theory for the algebras 7r and A'a to our permanent example — XXX magnet. 

Theorem 4.1 Let the representation T{u) of Tr be given by the product 

T{u) = KLn{u) . . .L2{u)Li{u) 

and the following nondegeneracy conditions be fulfilled. The first one concerns the matrix 
K 

det K 7^ K12 ^ 
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and ensures the Condition 3.1 of Lecture 3. The second condition coincides with the 
noninter section condition ( \3.1(^ for the sets {An}^^i defined by the formula l \3.1^ ) with 



— i 

Then all the conditions 3.1-4 of Lecture 3 are satisfied. The spectrum of operators 
{xn}n=i denoted now as X is the union of parallelepipeds 

X= U aX 

taken over all permutations a of coordinates {xn}^^i. The set X is defined by the 
formula i\3.1'li ). Note that the sets aX do not intersect due to condition l \3.1(^ ). The 



corresponding representation of the algebra X/^ is the direct sum of N\ standard irre- 
ducible representations with spec{xn}n=i = o"X described in the end of previous lecture, 
the parameters being defined above and ^± being arbitrary nonzero numbers satisfying 



The proof is performed by induction in A^. Let N = P + Q and 

Tiu)=Tpiu)TQiu) 

(A B\^(Ap Bp\(Aq Bq\ 
[C D ) [Cp Dp )[Cq Dq ) 

Let the generators of A'a corresponding to T,Tp,Tq be, respectively, zZ"^, xX"^, 
yY^ and the corresponding functions A±(n) be A±{u), A^^\ A|^\ 
Consider the eigenvalue problem 

for the matrix element B{u) of T{u) 

B{u) = Ap{u)Bq{u) + Bp{u)Dq{u) 

Note that since the spaces SymFunX ® SymFunY and FunX ® FunY are isomorphic 
one can think of $ as belonging to FunX ® FunY = Fun(X x Y). 

Substitutions u = Xp {p E {\, . . . , P}) and m = j/^ (g G {1, . . . , Q}) result, respec- 
tively, in 

/?(a;p)$(x,y) = 6q(xp - yi) . . . {xp - yQ)A^P{xp)^{E-^,y) 

and 

/3(l/,)$(x, y) = hp{y, - xi) . . . {y, - xp) (y,)$(x, E+y) 

27 



where 

p=i 



Let us extract out of $ the factor p 

*(x,y) = p(x,y)^^(x,y) 

satisfying the equations 

bgi^p -yi)--- i^P - yQ)piEp^, y) = p(x, y) 

bpiVq -xi)... {vg - xp)p(x, E+y) = p(x, y) 

It is easy to verify that the equations are compatible and have unique solution up to a 
constant coefficient. 

The original spectral problem is written in terms of as follows 

/3(x^)*(x,y) = A[^\x,)^{E;x,y) 

/5(|/,)vl/(x,y) = (|/,)vl/(x,i^+y) 
and apparently allows the separation of variables 

*(x,y)= {{[e.ix,)] (flxM 

\p=l J \q=l 

(5{x)ep{x) = A^^r^Opix -Tj) xe Af ) 
It remains to determine the spectrum of (3{u). Consider the equation 

(3{x)e{x) = A'L\x)e{x - 7]) 

There is the alternative: either ^(A~) ^ and hence (3{X~) = 0, or ^(A~) = 0. If 
0{\~) = then there is the new alternative: either 9{\~ +1]) ^ and then P{X~+rj) = 0, 
or 9{\~ + 1]) = 0, and so on. This argument shows that there exists A e A such that 
/9(A) = 0. Consequentely, 

(3{u) = (fliu - A,) j (fliu - \) j A, e Af ) A, e Af 
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or 

spec{zn}^^i = Z= [j aZ Z = XxY 



The spectrum of Zn being the same as expected, it remains to refer to the Theorem ^ 
claiming the uniqueness of the representations of the algebra corresponding to the 
sets crZ. 

Let us return now to our main problem — the spectral analysis of t{u) = A{u)+D{u). 
Consider the spectral problem 

t{u)ip = t{u)(p 

then substitute from the left u = Xn and use the definitions (|3.71 ), ( |3.7|) of together 
with the expression (|3.11|) for X^. The resulting set of equations 



TixrrM^) = A+(x)y;(E-x) + A;(x)^(E+x) {!,..., N} 

allows, obviously, separation of variables 

N 

ip{Xi, ...,Xn)= n QnM 



n=l 



which leads to the set of N one-dimensional finite-difference multiparameter spectral 
problems 

T{Xn)QniXn) = /^+iXn)QniXn-r]) + ^-{Xn)QniXn + 'n) ^^n ^ A„ U G {1, . . . , N} 

(4.1) 

which are identical in the form with Baxters's equation ( |3.3| ) arising in ABA. However, 
interpretation of the same equation is different. In the ABA case Q{u) was polyno- 
mial and now Q's are functions on discrete sets A„. In the ABA the natural way to 
solve the equation is to use the Bethe equations for the zeroes of Q{u). Though 
Bethe equations as a rule cannot be solved exactly they are simplified substantially 
in the infinite-volume limit N oo which is the most interesting case for applica- 
tions. In contrast, the FBA result suggests solving the system of recurrence relations 
for Q{u) numerically which can be of interest for small A^. So, for the XXX magnet 
ABA and FBA approaches are complementary. The FBA interpretation of the spectral 
problem ( |4.1| ) has however some advantage: since we have established the one-to-one 
correspondence between the eigenfunctions if of the original multidimensional spectral 
problem and those of the related one- dimensional ones there is no completeness prob- 
lem in FBA which arises in ABA due to the fact that, generally speaking, there could 
exist nonzero polynomial solutions to ( [4.1[ ) having no counterpart among eigenvectors 
of t{u) and vice versa there could be eigenvectors of t{u) which cannot be expressed 
as B{vi) . . . B{vm) |0). To be just, it is necessary to remark that the incompleteness of 
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ABA is a matter of degeneration and can always be removed by a small variation of 
parameters. 

Let us discuss now the possibility to apply FBA to the models other than XXX 
magnetic chains. A possible way to generalize the above results is to consider infinite- 
dimensional representations of the algebra Tr. If the spectrum X of the operators 
Xn is discrete then the generalization presents no difficulty. Consider for example the 
Goryachev-Chaplygin top discussed already in the previous lecture. It turns out that 
in order to construct the matrix T{u) representing the algebra it is necessary to 
introduce additional dynamical canonical variables p and q 

[p, q] = 1 

commuting with the generators Xa and Ja of e(3). Then the elements of T{u) are defined 
as follows 

A{u) = b{x+u- ^{x3,J+}) 

B{u) = e-^'^{u^-2J,u-{Jl + \)) 

C{u) = be''"^[{x+u-l{xs,J+}){u + p + 2J3)-bxl] 

D{u) = {u + p + 2J^){u^ -2J^u-{J^ - Jl + \) + hx-U-\{x^,J^}) 

where {, } is anticommutator, x± = Xi±ix2, J± = Ji±iJ2- The corresponding R-matrix 
is R{u) = u — 2V. It follows from the above formulas that 

t{u) = A{u) + D{u) = u^ + pu^ - 2{Hp + -)u - Gp 

8 

where 

Hp = H + pJs Gp = G + p{J^ - Jl + ^) 

H and G being the integrals of motion of the quantum GC top ( p.4|) , (|3l5|). Since p 
is one of integrals of motion it can be considered as a scalar parameter in Hp and Gp. 
In particular, for p = the GC top is recovered. The general case p 7^ has also nice 
physical interpretation: it corresponds to the so-called gyrostat. The zeroes of Mi,M2 of 
B{u) and the corresponding separation of variables are the same as discussed already 
in the previous lecture. 

Much more difficult is the case when the zeroes Xn of B{u) have continuous spectrum. 
Such an example is provided by the periodic Toda chain. In that case the spectrum of Xn 
is real and continuous but the shift t] in the Baxter's equations is purely imaginary. The 
rigorous mathematical justification of FBA presents in this situation serious analytical 
difficulties which are not yet overcome. However, quite formal application of FBA, "on 
the physical level of rigour" , leads to reasonable results and agrees with the results for 
the infinite volume limite obtained by independent methods. 

Another challenging problem is the generalization of FBA to Yangians of simple 
Lie algebras other then sl{2) that is those of rank > 1. I believe that solution of this 
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problem will clarify the algebraic roots of FBA scheme which at present is nothing that 
a misterious prescription: "take zeroes of B{u) and you obtain what you want". Let me 
conclude my lectures with a brief discussion of the problem for the sl{3) case. 

However, let us return first to the s/(2) case and look once more on the Baxter's 
equation ( p.3| ) dividing it by Q{x) 



t{u) = A+(uj — h A_(u 



Qiu) ' ' Q{u) 

Q{u±r^) 



Q{u) 



The formula can be rewritten as 

r(u) = A+(n) + A_(u) A±(n) = A^(u 

It is easy to see that 

AH = A_(n+|)A+(^-|) 

The fact that the quantities A±{u) are expressed in terms of the spectral invariants of the 
matrix T{u), its trace t{u) and quantum determinant A(u), suggests the interpretation 
of A-i-(m) as the eigenvalues of some operators which could be called "the quantum 
eigenvalues of the 2x2 matrix T(m)". One can write down also the quantum analogs 
of the secular equation for A±(m) 

A+(m)A+(u + r/) - t{u + ii)A+{u) + A(u + |) = 

A_(m - ri)A_{u) - t{u - r])A_{u) + A{u - |) = 
which turn into the familiar quadratic equation 

A^ - rA + A = 

in the classical limit rj ^ 0. Note that, in contrast with the classical case, A_|_ and A_ 
satisfy two different equations. 

The operators also can be considered as "quantum eigenvalues" of T(x„). Note 
that after subsitution u = Xn the matrix T{u) becomes triangular due to B{xn) = 0. In 
the classical case it means that the diagonal elements X~ = A{xn) and X+ = D{xn) 
of T{u) coincide with its eigenvalues. In the quantum case one can also write down the 
quantum analogs of the secular equations for X^. They are obtained by excluding X+ 
or X~ from the equalities 

t{Xn) = X+ + X„ 

A(a;„ ± |) = X^X^ 

and read 

[X+]2-X+t(a;J + A(x„ + |) = 
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The above observations need, of course, more profound study aimed at understanding 
the algebraic meaning of the "quantum eigenvalues" . However, they have already enough 
heuristic power to make some conclusions concerning FBA for sl{3) case. 

Consider the algebra 7r generated by the 3x3 matrix 



T{u) 



/ Tii(n) Ti2(u) Ti3(n) 
T2i(u) T22(m) T23(n) 
V Tsiiu) T32(m) T33(n) 



and the quadratic relations ( [L . 1|) with the R-matrix 

R(u) =u + r]V 

in ® C^. The algebra is equivalent to the Yangian 3^[s/(3)]. 
It turns out that the trace ti{u) of T{u) 

ti{u) = trT(u) 

does not provide the complete set of integrals of motion and should be accompanied 
with the operators 

t2{u)=tTi2Pi2T{u)f{u + r]) 

where P12 = {1 — V)/2 is the antisymmetrizer in ® C^, see (|1.4|) . 

The Casimir operator of Tr (quantum determinant of T{u)) is given by the formula 

A{u) = q-detT(M) = tri23Pr23T(M)f (u + ri)f{u + 2r]) 

where P]^3 is the antisymmetrizer in ®C^- 

The commuting quantities t2{u) and A('u) constitute three spectral invariants 

of the matrix T{u). 

The ABA for sl{N) case was developed in the papers by Yang, Sutherland, Kulish 
and Reshetikhin. It is the results of the last two authors which will be especially useful 
for our purposes. As Kulish and Reshetikhin (1982) have shown, the eigenvalues ti 2(m) 
of ti^2{u) together with the quantum determinant A(m) can be written down in the form 



Ti{u) = Ai(m) + A2(m) + A3(m) 

r2(n) = Ai{u)A2iu + r])+Ai{u)A3{u + T])+A2iu)As{u + r]) (4.2) 
A{u) = Ai{u)A2{u + ri)A3{u + 2ri) = Ai{u)A2{u + r])A3{u + 2ri) 

where the number polynomials Ai 2,3(m) like Aj-{u) for sl{2) case are expressed in terms 
of representation parameters. 
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The three "quantum eigenvalues" Ai^2,3(w) of T{u) can be expressed in terms of two 
polynomials Qi,2(w) 

Ai m =Ai M — — — A2 w = A2 w ^ / x TTT^ 

Qi{u) Qi[u) Q2[u) 

A3 M = A3 M — — — 
Q2{U) 

Kulish and Reshetikhin have derived also from ( |4.2| ) the finite-difference equations 
(analog of Baxter's equation) for Qi,2{u) which with the use of the shift /multiplication 
operators 

Si = Ai(x)e''^/^^ 
S2 = Ai(x-r7)A2(x)e''^/^^ 

can be put into the form 

[El-Ti{x + 2r])El + T2ix + 7])Ei- A{x)]Qi{x) = . . 

[El-T2{x + 7])El + Ti{x + 7])A{x)E2-A{x-r])A{x)]Q2{x) = ^ ' ' 

The last equations resemble very much the FBA separated equations and suggest, by 
analogy with the sl{2) case, the following conjecture concerning the possible form of 
FBA in the sl{3) case. 

Conjecture. There exist two sets of operators x^[^^\x[^^^ ni G {!,..., A^} and 
X2^^\X2^^^ ?T,2 G {1, . . . , 2A^} such that Xi^2 satisfy the "secular equations" 

X'l - X'^tiixi) + Xit2{xi) - A(xi) = Vni 

X| - X^t2{x2 -v)+ X2A{x2 - r])ti{x2) - A{x2 - r/)A(x2) = Vn2 

which lead to the separated equations ( ^T^) for the eigenvalue problem for the operators 
ti,2{u) ci'nd the eigenfunction 

N 2N 

^(xi,x2)= n n 

ni=l n.2=l 

The problem which remains unsolved is how to construct such operators that is to 
find a generalisation of the recipe B{u) = valid for the s/(2) case. In my opininon, 
solution of this problem will contribute to better understanding of the algebraic roots 
of Bethe Ansatz. 
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